Dysregulation and prometastatic function of glycosyltransferase C1GALT1 modulated by cHP1BP3/ miR-1-3p axis in bladder cancer

Background Abnormal glycosylation in a variety of cancer types is involved in tumor progression and chemoresistance. Glycosyltransferase C1GALT1, the key enzyme in conversion of Tn antigen to T antigen, is involved in both physiological and pathological conditions. However, the mechanisms of C1GALT1 in enhancing oncogenic phenotypes and its regulatory effects via non-coding RNA are unclear. Methods Abnormal expression of C1GALT1 and its products T antigen in human bladder cancer (BLCA) were evaluated with BLCA tissue, plasma samples and cell lines. Effects of C1GALT1 on migratory ability and proliferation were assessed in YTS-1 cells by transwell, CCK8 and colony formation assay in vitro and by mouse subcutaneous xenograft and trans-splenic metastasis models in vivo. Dysregulated circular RNAs (circRNAs) and microRNAs (miRNAs) were profiled in 3 pairs of bladder cancer tissues by RNA-seq. Effects of miR-1-3p and cHP1BP3 (circRNA derived from HP1BP3) on modulating C1GALT1 expression were investigated by target prediction program, correlation analysis and luciferase reporter assay. Functional roles of miR-1-3p and cHP1BP3 on migratory ability and proliferation in BLCA were also investigated by in vitro and in vivo experiments. Additionally, glycoproteomic analysis was employed to identify the target glycoproteins of C1GALT1. Results In this study, we demonstrated upregulation of C1GALT1 and its product T antigen in BLCA. C1GALT1 silencing suppressed migratory ability and proliferation of BLCA YTS-1 cells in vitro and in vivo. Subsets of circRNAs and miRNAs were dysregulated in BLCA tissues. miR-1-3p, which is reduced in BLCA tissues, inhibited transcription of C1GALT1 by binding directly to its 3′-untranslated region (3′-UTR). miR-1-3p overexpression resulted in decreased migratory ability and proliferation of YTS-1 cells. cHP1BP3 was upregulated in BLCA tissues, and served as an miR-1-3p “sponge”. cHP1BP3 was shown to modulate migratory ability, proliferation, and colony formation of YTS-1 cells, and displayed tumor-suppressing activity in BLCA. Target glycoproteins of C1GALT1, including integrins and MUC16, were identified. Conclusions This study reveals the pro-metastatic and proliferative function of upregulated glycosyltransferase C1GLAT1, and provides preliminary data on mechanisms underlying dysregulation of C1GALT1 via miR-1-3p / cHP1BP3 axis in BLCA. Supplementary Information The online version contains supplementary material available at 10.1186/s13046-022-02438-7.


Background
Glycosylation is a complex and commonly occurring type of protein modification that promotes functional plasticity by increasing protein heterogeneity [1,2]. Normal O-glycosylation, one form of glycosylation, is essential for proper development [2] and cellular differentiation [3]; on the other hand, dysregulation of O-glycans contributes to a variety of pathogenic processes, including inflammation [4], tumorigenesis [5], and metastasis [6]. Addition of N-acetylgalactosamine (GalNAc) to serine or threonine, catalyzed by polypeptide GalNAc transferases (ppGalNAcTs), results in formation of Tn antigen (GalNAcα-O-Ser/Thr). O-glycans based on Tn antigen are extended to various complex, branched structures by sequential glycosyltransferase reactions. Core 1 β1,3-galactosyltransferase (C1GALT1; also termed T-synthase) transfers galactose (Gal) to GalNAc to form Galβ1-3GalNAcα1-Ser/Thr (core 1 O-glycan structure, termed T antigen). Core 1 structure is the precursor for subsequent extension and maturation of mucin-type O-glycans.
Expression of truncated O-glycans, such as sialyl-Tn (STn) and Tn antigens, is frequently observed in gastric, hepatic, and some other types of cancer, and is correlated with cancer aggressiveness and progression [7,8]. Truncated O-glycans alter antigenic characters and physiological functions of glycoproteins, and promote invasiveness, chemotherapy resistance and metastasis of cancer cells [9,10]. Some of them have been developed as small-molecule inhibitors [11] and carbohydrate vaccines for application in cancer therapy [12]. Better understanding of how truncated O-glycans are regulated may therefore lead to improved immunotherapeutic strategies [13].
C1GALT1 is localized in Golgi apparatus, and its enzymatic activity depends on an unusual molecular chaperone termed "Cosmc" (core 1 β3GalT specific molecular chaperone). C1GALT1 is commonly upregulated during tumorigenesis, including development of hepatocellular carcinoma [8], breast cancer [14], and some other types of cancer [15]. C1GALT1 has a modulatory effect on expression of target glycoproteins such as MUC16 [16], Muc5Ac, and MUC1 [17], and regulates progression of certain types of cancer. In gastric cancer, it promoted phosphorylation of EPHA2 and enhanced soluble Ephrin A1-mediated migration by modifying EPHA2 O-glycosylation [18]. In pancreatic adenocarcinoma, its knockdown suppressed focal adhesion kinase (FAK) phosphorylation at Y397/Y925, and reduced cell/ extracellular matrix adhesion by modifying integrin αv [15]. A loss of C1GALT1 resulted in a Tn enrichment on CD44, which activated ERK/NF-κB signaling and facilitated the cancer stem cell features in pancreatic cancer cells [19]. However, C1GALT1 functioned differentially in other studies. It is more highly expressed in well-differentiated than in poorly-differentiated pancreatic cancer tissues [16]. Modulation of C1GALT1 expression contributes to alterations of tumorigenesis and metastasis. Mice with knockdown of gastric epithelial C1GALT1 developed spontaneous chronic gastritis and gastric cancer [17]. Downregulation of C1GALT1 suppressed the TrkA expression and promoted malignant behaviors of neuroblastoma cells [20]. To date, few studies have addressed the role of C1GALT1 in bladder cancer (BLCA), and molecular mechanisms that underlie C1GALT1 dysregulation are unclear.
The present study revealed abnormal C1GALT1 levels in BLCA. Normal bladder tissues, in comparison with BLCA tissues, showed differential expression of micro-RNAs (miRNAs) and circular RNAs (circRNAs). We provide preliminary data on mechanisms underlying dysregulation of C1GALT1 via mRNA/ miRNA/ circRNA axis, and effects of C1GALT1 and corresponding regulatory non-coding RNAs on cell proliferation and migratory ability.
All cells were used for the experiments between passages 3 and 15.

Tissue microarray (TMA) analysis
TMA analyses were performed as described previously [22]. TMAs containing BLCA and paracarcinoma tissues were from Shanghai Outdo Biotech Co. In brief, TMAs were deparaffinized, rehydrated, subjected to antigen retrieval, blocked, incubated with biotinylated lectin or primary antibody overnight at 4 °C, rinsed with 1× PBS (0.01 mol/L phosphate buffer containing 0.15 mol/L NaCl, pH 7.4), incubated with HRP-streptavidin or secondary antibody for 1 h at 37 °C, visualized with DAB reagent, stained with hematoxylin, and photographed. Mean optical densities of the TMA staining signals were calculated using Image-Pro Plus software (Media Cybernetics; Silver Spring, MD, USA).

Knockdown of C1GALT1 and cHP1BP3
shRNA targeting C1GALT1 or backsplice junction of cHP1BP3 were respectively subcloned into lentivirus vector (Tet-pLKO-puro vector, Addgene; Cambridge, MA, USA) to construct sh-C1GALT1 and sh-cHP1BP3 vector, and shRNA sequences were listed in Table S1. Constructed lentiviral vector and two assistant vectors (pMD2.G and psPAX2) were transfected into HEK293T cells with PEI MAX 40000 (Polysciences; Warrington, PA, USA). Virus particles were collected after 2 days, applied to infect YTS-1 cells, and stable transfectants were selected using puromycin. Transient silencing of C1GALT1 and cHP1BP3 in T24 was performed by transfection of constructed vectors. Silencing of C1GALT1 and cHP1BP3 was induced by doxycycline (Doxy).

Proliferation assay
Cell proliferation was assayed as described previously [25]. For Edu assay, cells treated with/ without Doxy were incubated with 50 μM Edu for 6 h, fixed with 4% paraformaldehyde, permeabilized with 0.5% Triton X-100, stained with iClick EdU solution (GeneCopoeia; Rockville, MD, USA), rinsed with 1× PBS, and subjected to flow cytometry. For CCK8 assay, cells in 96-well plates were incubated with CCK8 solution (Beyotime) for 4 h, and absorbance at 450 nm was determined using a microplate reader.

Colony formation assay
Colony formation was assayed as described previously [25]. In brief, cells in 6-well plates (1000 cells/ well) were cultured for 1-2 wk. until colonies were clearly observable. Colonies were rinsed with 1× PBS, fixed with 4% paraformaldehyde, stained with crystal violet solution, and photographed. Acetic acid solution was added to dissolve crystal violet, and absorbance at 595 nm was determined.

Transwell assay
This assay was performed as described previously [22]. 10 4 cells suspended in FBS-free medium were inoculated into upper chamber, and medium supplemented with 10% FBS was added into lower chamber. After 24 h culture, cell migrated across the membrane were fixed, stained with crystal violet, and photographed under microscope.

Patient samples
Sera and tissues from normal subjects and BLCA patients were obtained from the Third Affiliated Hospital of Xi'an Jiaotong University. Written informed consent was obtained from all patients, in accordance with the Declaration of Helsinki. Experiments using human tissues were approved by the Research Ethics Committee of Northwest University.

RNA sequencing analyses
For circRNA sequencing analysis, total RNAs were isolated from BLCA and adjacent normal tissues using TRIzol reagent kit (Invitrogen) as per manufacturer's instructions, treated with RNase R, and purified using RNeasy MinElute Cleanup Kit (Qiagen; Valencia, CA, USA). Strand-specific library was constructed using VAHTS Total RNA-seq Library Prep Kit for Illumina, as per manufacturer's instructions. In brief, circRNAs were enriched by removal of rRNAs, fragmented into short fragments, and reverse-transcribed into cDNA. cDNA was further amplified, and sequenced using Illumina HiSeq 2500 by Gene Denovo Biotechnology (Guangzhou).
For miRNA sequencing analysis, isolated total RNAs (size range 18-30 nt) were enriched by PAGE. 3′ and 5′ adapters were ligated to RNAs, and ligation products were reverse-transcribed by PCR amplification, enriched to generate a cDNA library, and sequenced as above.

RNase R treatment
RNAs from YTS-1 cells were treated with 3 U/μg RNase R (Geneseed; Guangzhou), incubated for 30 min at 37 °C, reverse-transcribed with primers, and detected by reverse transcription quantitative real-time PCR (RT-qPCR), and primers were listed in Table S1.

Identification of proteins with T antigen
Proteins were denatured, reduced, alkylated, digested, and desalted as described previously [22]. Lyophilized peptides (150 μg) were dissolved with 50 mM NH 4 HCO 3 , and incubated with 75 μL PNA-agarose (Vector Labs) overnight at 4 °C. The mixture was rinsed with PBS, and boiled for 10 min to release peptides. Glycopeptides with T antigen were collected by centrifugation and purified using Oasis HLB cartridges (Waters; Milford, MA, USA). Two-dimensional LC-MS/MS and data analysis were performed using LTQ Orbitrap MS (Thermo Fisher), and MaxQuant software program with T antigen as variable modification, as described previously [29,30].

Tumor formation in mice
Animal experiments were performed in accordance with guidelines of the Animal Care and Use Committee of Northwest University. Effects of C1GALT1 on BLCA tumor proliferation and metastasis were studied using mouse subcutaneous (s.c.) xenograft, trans-splenic metastasis and popliteal lymphatic nodes (LNs) models.
The role of C1GALT1 in BLCA based on blood metastasis to liver was investigated by injecting BLCA cells (10 5 cells in 100 μL PBS) into spleens of 4-to 6-wk-old mice as described previously [31]. After 4 wk., mice were euthanized, and tumors were excised and weighed.
Tumors, spleens, and livers from each type of experiment were dissected, fixed, and paraffin-embedded for histopathological analysis.
PLNs metastasis model was performed by the injection of 100 μL PBS suspensions of BLCA cells (10 5 cells in 100 μL PBS) into footpads of 4-to 6-wk-old mice as described previously [32]. After 4 wk., popliteal lymphatic nodes were dissected, fixed, and paraffin-embedded for histopathological analysis.

Patient-derived xenografts (PDXs) models
PDX tumor was derived from a male patient with papillary urothelial carcinoma. Briefly, three serial passage xenograft tumor pieces of ~ 60 mm were subcutaneously grafted into the 4-to 6-wk-old mice. When tumors grown to ~ 200 mm, the mice were randomly divided into vehicle or itraconazole (ITZ) subgroups. Solvent control or ITZ (20 mg/mL, twice daily) was given orally. Tumor size was measured, after which tumors were excised, weighed and subjected to histopathological analysis.

Statistical analysis
Data are presented as mean ± SD from three independent experiments unless specified otherwise. Two-tailed Student's t-test was used for comparison of data sets between two groups, and differences with p < 0.05 were considered statistically significant. Notations in figures: *p < 0.05; **p < 0.01; ***p < 0.001.

Upregulation of C1GALT1 expression in BLCA
Expression of C1GALT1 and its product T antigen (Fig. 1A) were evaluated by TMA analysis in 44 pairs of primary human BLCA tissue samples with matched adjacent noncancerous bladder tissues (Table S2). C1GALT1 expression was clearly higher in BLCA than in noncancerous cells (Fig. 1B), as were corresponding T antigen levels (Fig. 1C). Western and lectin blotting likewise revealed increased C1GALT1 expression and T antigen levels in BLCA tissues (Fig. 1D), and higher serum levels of T antigen in BLCA tissues were shown by lectin PNA-based ELISA (Fig. 1E). The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases indicate significantly higher C1GALT1 expression at mRNA level in BLCA relative to adjacent tissues (Fig. 1F, G, S1A). According to Kaplan-Meier curve analyses from TCGA database, overall survival prognosis was significantly worse for patients with high (vs. low) C1GALT1 expression (Fig. 1H), and high C1GALT1 expression in BLCA patients was correlated with poor disease-free survival (Fig. 1I). Similarly, high (vs. low) expression C1GALT1 or T antigen showed worse overall survival using TMA analyses (S1B, C).
In the present study, C1GALT1 expression was higher in invasive BLCA YTS-1 than in normal bladder HCV-29 cells (Fig. S1D). And O-glycan analysis showed higher T antigen levels in bladder cancer cells KK47, T24 and YTS-1, than in normal epithelial bladder cells HUC-1 and HCV-29  (Fig. 1J). These findings indicate that dysregulated expression of C1GALT1 and T antigen is characteristic of BLCA.
Tumor volume and weight in vivo were significantly reduced by ITZ treatment or C1GALT1 silencing (Fig. 2G-I). Immunohistochemical analyses revealed decreased T antigen and C1GALT1 expression and cell proliferation, and enhanced apoptosis, in ITZ-treated and C1GALT1-silenced tumor xenografts (Fig. S4A). We further established subcutaneous BLCA PDXs in mice. The tumors derived from PDXs with ITZ treatment regrew shortly compared to control group (Fig. 2J-L). Using PDXs model, we found that ITZ treatment inhibited C1GALT1 expression, enhanced apoptosis, and reduced proliferation (Fig. S4B).
In the trans-splenic metastasis model, C1GALT1 silencing significantly suppressed tumor metastasis to liver (Fig. 2M, N), enhanced apoptosis, and reduced proliferation (Fig. S5). In the popliteal lymphatic metastasis model (Fig. 2O), the volumes of lymphatic nodes were smaller in the C1GALT1 silencing group than which in control group (Fig. 2P). C1GALT1 silence suppressed the ability of the BLCA cells to metastasize to the LNs, as determined by the number of metastatic LNs and Haematoxylin & Eosin (H&E) staining (Fig. 2Q, S6A, B).
In vivo studies revealed significant reduction of tumor volume and weight in miR-1-3p-overexpressing cells (Fig. 4H-J). Immunohistochemical analysis of (See figure on next page.) these cells showed decreased C1GALT1 and T antigen expression and cell proliferation, and increased apoptosis (Fig. S7H). In the trans-splenic metastasis model, miR-1-3p overexpression significantly reduced tumor metastasis to liver and cell proliferation, and increased apoptosis (Fig. 4K, L, S8).

Screening regulatory circRNAs of miR-1-3p
CircRNAs have been shown to act as miRNA "sponges", and to play essential roles at the transcriptional level.
Association of circRNAs with miR-1-3p was investigated by characterizing circRNA transcripts through RNA-seq analysis of rRNA-depleted and RNase-treated RNA from pairs of BLCA tissues as described above. We identified 10,055 circRNAs in these tissues, including 4129 that overlapped with circBank [35] and 5926 novel circRNAs (Fig. 5A). Of these, 8438 circRNAs (83.9%) consisted of protein coding exons from 3519 host genes, and the rest consisted of exon-intron, intron, intergenic, or antisense (Fig. 5B). Most of the exonic circRNAs had lengths < 1000 nt (Fig. 5C). Identified circRNAs were quantified, and 414 were found to be differentially expressed in BLCA tissues relative to matched adjacent normal tissues (Table S4). Of these differentially expressed circR-NAs (termed DEcircRs), 283 were upregulated and 131 were downregulated in BLCA tissues (Fig. 5D). PCA of DEcircRs demonstrated clear distinction of those in BLCA vs. adjacent normal tissues (Fig. 5E). Of the 414 DEcircRs, 349 were exonic, of which 244 were documented in circBank. For confirmation of RNA-seq findings, we examined the top four upregulated and bottom three downregulated DEcircRs. Head-to-tail splicing was assayed by RT-qPCR, with divergent primers and Sanger sequencing (Figs. 5F-G, S9A-B). Head-to-tail junctions of the seven DEcircRs as above were all confirmed, demonstrating high accuracy of our circRNA profiling. We examined sensitivity of the seven DEcircRs to RNase R digestion, to rule out possible production of head-to-tail splicing by trans-splicing, genomic rearrangements, or PCR artefacts. The DEcircRs all showed greater resistance to RNase R than that of GAPDH control (Fig. 5H).
In conclusion, characterized circRNAs were consisted mainly of exons, and seven selected DEcircRs were confirmed as having head-to-tail junctions.

Sponage role of cHP1BP3 on miR-1-3p
Potential regulatory relationships between miR-1-3p and circRNAs were evaluated using the three public databases. 213 circRNAs were predicted to target hsa-miR-1-3p based on our results (Fig. 6A). In view of the miRNA sponge function of circRNA, 14 circRNAs were filtered by overlap between DEcircRs and predicted target circRNAs of miR-1-3p (Fig. 6B). Among these, the top four upregulated circRNAs were selected for examining interaction with miR-1-3p. Biotin-labeled miR-1-3p was introduced into YTS-1 cells, and interactions were evaluated by RNA pulldown, which revealed strong binding of circRNA cHP1BP3 to miR-1-3p (Fig. 6C). Luciferase reporter assay indicated direct binding of miR-1-3p to predicted cHP1BP3 binding sites (Fig. 6D), other than mutated sites (Figs. 6E, S10A). cHP1BP3 was highly stable in response to transcriptional blockage, other than linear counterpart GAPDH (Fig. 6F). cHP1BP3 was significantly downregulated relative to oligo (dT) 18 primers when random hexamer primers were used for reverse transcription experiments (Fig. 6G), indicating the circular nature of cHP1BP3. FISH assay showed that cHP1BP3 was located mainly in cytoplasm, and significantly enhanced in BLCA vs. adjacent normal tissues (Fig. S10B). cHP1BP3 level was negatively correlated with linear counterpart (Fig. S10C), indicating that its higher expression in BLCA reflected its functionality, and was not simply a by-product of splicing. Moreover, cHP1BP3 expression was significantly increased in BLCA tissues compared to normal tissues (Fig. 6H). cHP1BP3 was silenced in YTS-1 and T24 cells to further elucidate its physiological functions. cHP1BP3 expression was significantly reduced by Doxy treatment using divergent primers (Fig. S11A), and cHP1BP3 silencing suppressed T antigen and C1GALT1 expression in YTS-1 and T24 cells (Fig. S11B, C). cHP1BP3 silencing also reduced cell proliferation, migratory ability, colony formation, and doxorubicin resistance relative to vehicle control (Figs. 6I-K, S11D-F), and suppressed tumor volume and weight in vivo (Fig. 6L-N). Immunohistochemical analysis confirmed that cHP1BP3 silencing suppressed T antigen and C1GALT1 expression, leading to reduced cell proliferation and increased apoptosis (Fig.  S11G). In the trans-splenic metastasis model, cHP1BP3 silencing significantly suppressed tumor metastasis to liver (Fig. 6O, P), enhanced apoptosis, and reduced proliferation (Fig. S12), similarly to miR-1-3p overexpression (Fig. 4).

Target proteins of C1GALT1
C1GALT1 controls the essential step of O-glycosylation elongation, and strongly affects target protein function. Target O-glycosylated proteins of C1GALT1 involved in BLCA were examined by glycoproteomic analysis in combination with PNA enrichment and triplicate LC-MS/MS (Fig. 7A). 117 glycopeptides were identified in C1GALT1silenced and parental YTS-1 cells, covering 67 glycoproteins (Table S5). PCA showed clear distinction of these (See figure on next page.) Fig. 6 Effects of cHP1BP3 on malignant behavior of BLCA cells. A Venn diagram of predicted circRNAs acting as miR-1-3p sponges, from miRanda, RNAhybrid, and TargetScan databases. B Venn diagram of identified upregulated circRNAs and predicted circRNAs acting as miR-1-3p sponges. C RT-qPCR analysis of circRNA levels in streptavidin-captured fractions from YTS-1 cell lysates following transfection with 3′-end biotinylated miR-1-3p or control RNA (NC). **, p < 0.01. D Confirmation by luciferase reporter assay of miR-1-3p-binding sites on circRNAs. Control miRNA or miR-1-3p mimics were co-transfected with luciferase reporters containing predicted miR-1-3p binding sites. *, p < 0.05; **, p < 0.01; ***, p < 0.001. E Luciferase activity of luciferase-cHP1BP3 or its mutant in HEK293 cells following co-transfection with miR-1-3p. ***, p < 0.001. F Stability of cHP1BP3 following glycoproteins in C1GALT1-silenced vs. parental YTS-1 (Fig. 7B). Gene Ontology (GO) and KEGG enrichment analyses revealed enrichment of proteins involved in protein transport from Golgi to plasma membrane, Rho GTPase cycle, and vesicle localization (Fig. 7C). Proteinprotein interaction network analysis indicated that target proteins formed clusters with enriched histone kinase activity and calmodulin binding (Fig. 7D). Complete hierarchical clustering and visualization (heatmap) showed clustering of the identified glycoproteins (Fig. 7E), and their expression patterns in C1GALT1-silenced and parental YTS-1 cells. Quantification of the glycopeptides showed that three (including MUC16 and HMCN1) were differentially expressed (Fig. 7E). Consistently, immunoprecipitation analysis indicated that MUC16 was decorated with T antigen and MUC16 expression was decreased in C1GALT1 silenced YTS-1 cells (Fig. 7F). C1GALT1 silencing facilitated the degradation of MUC16 in YTS-1 cells treated with cycloheximide, the inhibitor of protein synthesis (Fig. 7G). MUC16 silencing by siRNA (Fig. 7H) suppressed the cell proliferation and migratory ability of YTS-1 cells (Fig. 7I, J). Collectively, these results indicated that C1GALT1 silencing suppressed migratory and proliferation of BLCA cells by modifying target glycoproteins including MUC16.

Discussion
Truncated O-glycans and responsible enzyme C1GALT1 play important roles in a variety of physiological and pathological processes, including cancer development and progression [36,37]. C1GALT1 overexpression has been reported for various cancers [18,[38][39][40]. Recent studies demonstrated that T antigen, the product of C1GALT1 activity, was present in high-grade BLCA but not in normal urothelium [41]. BLCA cells under hypoxic conditions showed striking upregulation of C1GALT1 and associated downregulation of C2GnT, the enzyme that catalyzes O-glycan core 2 formation with T antigen as substrate; this observation may account for T antigen accumulation in BLCA [41]. Consistant with other reports, we found that T antigen and its responsible enzyme C1GALT1 were up-regulated in BLCA, and C1GALT1 expression was positively correlated with malignant phenotypes by modulating C1GALT1 expression in BLCA cells. Notably, the treatment of ITZ, as one antifungal medication, could down-regulate the C1GALT1 expression, and suppress BLCA tumor growth and migratory ability in cell-and patient-derived xenografts, indicating the therapeutic potential of ITZ in BLCA.
In recent studies, C1GALT1 has been shown to be negatively modulated, similarly, by miR-181d-5p in lung adenocarcinoma (19), miR-152 in gastric cancer (16), and miR-124-3p in aging colon (20). In this study, C1GALT1 was revealed to be regulated by miR-1-3p, not by other miRNAs. These findings indicate that miRNAs targeting C1GALT1 have differential functions depending on cancer type. In a study of prostate cancer, miR-1-3p showed downregulation, and its ectopic expression suppressed cell growth and cell cycle progression in vitro and in vivo through targeting of 3′-UTRs of two key cell cycle genes, E2F5 and PFTK1 [42]. In non-small cell lung cancer, levels of miR-1 were low, and its overexpression inhibited tumor growth and angiogenesis by targeting Mpl gene [43], similarly to our findings.
circRNAs have been shown to act as miRNA sponges, participate in target gene splicing, transcribe genes into proteins, and interact with RNA-binding proteins. We therefore examined the regulatory role of circRNA in abnormal C1GALT1 expression in BLCA. cHP1BP3 comprising several exons was identified as a significantly upregulated circRNA in BLCA, and was predicted by bioinformatic analysis to bind miR-1-3p. cHP1BP3 functions as a competing endogenous RNA (ceRNA), and was shown by biotinylated miRNA pulldown and luciferase reporter assays to bind competitively to miR-1-3p. In rescue experiments, cHP1BP3 abolished the endogenous suppressive effect of miR-1-3p on target C1GALT1. These observations, in combination with our identification of C1GALT1 as the direct target of miR-1-3p, indicate that C1GALT1 is regulated via cHP1BP3/ miR-1-3p axis in BLCA (Fig. 8).
C1GALT1 is able to modify O-glycans on target proteins and regulate their physiological functions. We performed PNA enrichment of T antigen-bearing glycoproteins followed by MS analysis, and observed that integrin α7 was bearing T antigen. Previous studies have demonstrated similar O-glycosylation changes on integrin β1, αv, and α5 subunits. C1GALT1 affects cell migratory ability and invasiveness by modifying O-glycans on integrin α5 [44] or β1 [39], and modulating their downstream signaling. C1GALT1 knockdown enhanced radiosensitivity of esophageal cancer cells, and suppressed irradiation-enhanced migration and invasion by modifying glycosylation of integrin β1 [45]. MUC16 has also been identified as a C1GALT1 target glycoprotein in BLCA. Pancreas-specific C1GALT1 knockdown resulted in early onset of tumors and metastatic lesions, possibly because of dysregulated MUC16 O-glycan modification [16]. Deletion of MUC16 decreases tumor cell migration in pancreatic ductal adenocarcinoma cells [46], which is in line with the decrease in proliferation and migratory ability in MUC16silenced YTS-1 cells. MUC16 bearing truncated O-glycans exacerbates malignancy by FAK activation through interactions with integrins on cell membranes [46], indicating that O-glycans might regulate the localization and turnover of target proteins and be involved in cancer progression.

Conclusions
In conclusion, our findings indicate that C1GALT1 is modulated by cHP1BP3/ miR-1-3p axis and upregulated in BLCA. C1GALT1 silencing suppressed migratory ability and proliferation of BLCA cells in vitro and in vivo by modifying target glycoproteins. Notably, C1GALT1 is revealed to be regulated via circRNA/ miRNA axis in BLCA, which might provide theoretical basis for the diagnosis and therapy of BLCA.